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Abstract 

We develop a supersymmetric representation of the Hubbard operator which 
unifies the slave boson and slave fermion representations into a single U{\) x 
SU gauge theory, a group with larger symmetry than the product of two 
U{1) gauge groups. These representations of the Hubbard operator can be 
used to incorporate strong Hund's interactions in multi-electron atoms as a 
constraint. We show how this method can be combined with the SP(N) group 
to yield a locally supersymmetric large-N formulation of the t — J model. 
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One of the fascinating aspects of strongly correlated materials is their propensity to de- 
velop novel metallic phases in situations where local moments interact strongly with mobile 
electrons. Examples of such situations include metals near a metal insulator transition,!! met- 
als at an anti-ferromagnetic quantum critical point)! and anti-ferromagnetic heavy fermion 
superconductors.! These discoveries challenge our understanding of how spin and charge 
interact at the brink of magnetism. 

Theoretical approaches to these problems are hindered by the difficulty of capturing 
the profound transformation in spin correlations that develops at the boundary between 
antiferromagnetism and paramagnetism. Usually we model these features by representing 
the spin as a boson in a magnetic phase,! or as a fermion in a paramagnetic phase,! but 
by making this choice, the character of spin and charge excitations which appear in an 
approximate field theory is restricted and lacks the flexibility to describe the co-existence of 
strong magnetic correlations within a paramagnetic phase. 

These considerations have motivated the development of new methods to describe the 
spin and charge excitations of a strongly correlated material which avoid making the choice 
between a bosonic or fermionic spin.iHl This paper attempts to stimulate further progress in 
this direction by introducing a supersymmetric representation of Hubbard operators. The 
method used here is an extension of the supersymmetric spin representation introduced by 
Coleman, Pepin and Tsvelikffi (CPT). Remarkably, the supersymmetry in the CPT spin 
representation survives the introduction of charge degrees of freedom, opening the method 
to a wider range of models. 

Hubbard operators^ provide a way to describe atoms in which Coulomb repulsion pre- 
vents double-occupancy of a given orbital. Suppose \a) G {|0), \a)} describes a set of 
atomic states involving a charged "hole" |0) or a neutral spin state |cr) with spin com- 
ponent a G {1 . . . N} which for generality can have one of N possible values. The Hubbard 
operators are written 

X ab =\a){b\ (1) 
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where a, b G {0,iV}, represent an atomic state with N possible spin configurations. The 
operators X aa i are bosonic spin operators whereas the X a0 and X 0a are fermionic operators 
that respectively create and annihilate a single electron. The spin operators X au i are the 
generators of the group SU(N). The additional operators expand the group to a supergroup 
SU (iV|l)0 that describes the physical spin and charge degrees of freedom of the atom. These 
operators satisfy a graded Lie algebra 

[X a b, X cd ]± = 5 bc X ad ± 5 a dX cb . (2) 

where the plus sign is only used for fermionic operators. The absence of a Wick's theorem 
for these operators is normally overcome by factorizing the fermionic Hubbard operators 
as a product of canonical creation and annihilation operators. This can be done by repre- 
senting the empty state by a "slave boson" and the spin by a fermioni or alternatively, by 
representing the empty state as a "slave fermion" and the spin by a Schwinger boson.0 
We now generalize this approach, introducing 

F a — (fi, . . . /at, (p) 

B a = (h,...b N , X ) (3) 

where b a and f a denote a Schwinger boson! and Abrikosov pseudo-fermion0 respectively, 
while and x are slave bosons! and fermions0 respectively. In terms of these operators, 
the supersymmetric representation of the Hubbard operators is written 

X ab = B\B b + FlF b , (4) 

Written out explicitly, this is 

X aa , = b%, + /t/ CT , 

X*o = blx + fa, X 0a = X % + tff* 

*oo=X f X + 0V (5) 

By summing the slave fermion and slave boson representations we are guaranteed that the 
representation satisfies the correct commutation algebra. The novelty of our approach lies 



in the two unique constraints which make the representation irreducible, which we show to 
be 



Q = n b + ri0 + n f + n x , (6) 

the total number of particles and 

Y = n 4> + n f -(n b + n x ) + -[9,9\ (7) 

the "asymmetry" of the representation, where 9 = J2a tfrfo ~ xV an d its conjugate 6^ are 
fermionic operators which satisfy the algebra {9, 9^} = Q. The 9 operators interconvert 
bosons and fermions. 

9^ 9^ 

fra^=^ fa, ~X^r <f> 

9 9 

The special feature of this representation is that 9 and 9^ commute with the constraints 
= [6"\ Y] = 0, the bosonic Hubbard operators 

[9«\X ca ,\ = [6®,X O0 ] = 0, 

and they also anti-commute with the fermionic Hubbard operators 

{9^,X a0 } = {9^,X 0a } = 0, 

so that there is a local supersymmetry which underlies the constraint. The operators Q, 9 
and 0t are the generators of the simplest supergroup the operator Y generates 

an additional £7(1) symmetry. Remarkably, by combining the slave boson and slave fermion 
representations, the abelian gauge groups of the starting representation "fuse" into a su- 
pergroup with greater symmetry Usb(1) x Usf(1) U(l) x S77(l|l). If we introduce the 

operator A = [f}9 — 9'rj\, where rj and fj are Grassman numbers, then under an S77(l|l) 

(B a \ 

rotation, the fields ip a — \ transform as 

\Fa) 

i> a - e A *p a e- A = ,/> a + [AM + ±[A,[A, A}] (8) 



where the Grassman coefficients truncate the expansion at second-order. Expanding this 
expression gives tp a — ► hi/) a , ip\ — > ijr h h) where 

/ y/l-fjfq -fj \ 

V 77 V 1 - vvj 

is an S77(l|l) matrix, satisfying h'h = 1 The X-operators (Q) can be written as X a b = iplipb- 
Under the action of the SU group, X ab — > iplh'hijjb = X a &, explicity demonstrating the 
local gauge invariance. 

To guarantee that the Hubbard operator representation is irreducible, we need to set 
the values of the linear and quadratic Casimirs of the SU(N\1) group. Under the SU(N\1) 
group, the spinors B and F transform according to B — > BU, F — > where U = U st 

denotes the supertranspose of the unitary SU(N\1) matrix, U 0.The Hubbard operators 
X ab = B\B h + F^F a thus transform according to X ab — > {U^XU) ah . Since U^U = 1, it 
follows that U st (W) st = 1. However, the supertranspose and hermitian conjugate do not 
commute and are related by (f/^) s * = g(U st Yg, where g = Diag[l ... 1, —1] is the invariant 
metric tensor of SU(N\1). Thus the U are not unitary, but satisfy UgU^ = g. Using the 
property that Tr[A5] = Tr[BgAg], it follows that 

= Tr[X], C {2) = Tr[XgX\, (9) 

are invariant under the transformation X — > WXU. These are the linear and quadratic 

Casimirs of the SU(N\1) group. Inserting (^) into (|9|), we find that C^ 1 -* = Q, while the 
quadratic Casimir is 

C*( 2 ) = X aa iX a t a — Xo-oXqo- + X 0r7 'X a ' — Xq (10) 

where summation over a, a' G {1, N} is implied. When we expand the Casimir in terms of 
the canonical creation and annihilation operators, we find that 

C {2) =Q(N -1-Y), (11) 

with Q and Y as given in (g) and (0 ). So by defining Y and Q, we uniquely set the 
representation. 
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FIG. 1. (a) Fundamental representation (Q,Y) = (1,0), (b) L-shaped Young tableau corre- 
sponding to the spin representation generated by supersymmetric Hubbard operators. The asym- 
metry Y = h—w and Q is the number of boxes, (c) Young tableau for fully symmetric representation 
corresponding to Slave fermion limit (d) Fully antisymmetric, slave boson limit. 

Each conserved value of (Q,Y) describes an irreducible representation of the SU(N\1) 
group; the fundamental representation, (Q, Y) = (1, 0) corresponds to an atomic orbital 
with no double occupancy (Fig. |l|(a)). More general representations involve spin wavefunc- 
tions with symmetric and antisymmetric correlations, denoted by an "L— shaped" Young 
tableau^! with Q boxes, where Y = h — w is the difference between the height and width 
(Fig. [I] (b-d)). These representations describe the physics of multi-electron atoms where 
the spins are Hund's coupled, and in this way strong Hund's couplings can be incorporated 
into an infinite U Anderson model using the constraints (||) and (]?]). As an example, the 
material LiX^O^ develops a paramagnetic heavy fermion ground-stateE! in which vanadium 
ions form a mixed valence admixture of a ^(S = 1/2) and a Hund's coupled d 2 (S = 1) state. 
Since the electrons in the d 2 configuration are in a symmetric wavefunction, corresponding 
to a row-tableau, this situation is described by Hubbard operators in the representation 
(Q,Y) = (2,-1): 
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As a second example, consider U Pd^Al-^ in which uranium atoms fluctuate between an f 2 and 
an f 3 configuration. Surprisingly, part of the spin magnetically orders, while the remainder 
forms a singlet superconductor with the conduction electrons.! In this case, the f-electrons 
are spin-orbit coupled, with j — 5/2, forming an SU(N) multiplet with N = 2j + 1 = 6. In 
practice, crystal field effects break this large degeneracy, but a toy model for the physics can 
be obtained using SU (N) Hubbard operators to describe the charge fluctuations, subject to 
the constraint (Q,Y) = (2, 1). This leads to valence fluctuations involving an L— shaped 
spin f 3 spin configuration: 

P f 3 

e~ + 



In this scheme the vertical leg of the representation can form a singlet with conduction 
electrons, leaving a single residual spin free to magnetically order.0 

In many problems we are interested in interacting atoms containing either one, or zero 
electrons. Physical states corresponding to this situation have Q = 1, Y = 0: 

Q|-0> = |-0>, Y\^) = 0. (12) 

These conditions do not force the representation into a slave boson, or slave fermion repre- 
sentation. Here, it is useful to note that 9 and 9^ behave as lowering, and raising operators. 
In fact, because {9, 9^} = Q, 

1 .4 1 - 1 
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T- = -=6, r z = [r + ,r_] = -[9\9], 



VQ VQ 2 L+ ' J Q 

behave as the raising, lowering and z components of a "superspin" operator. If we take the 
sum and difference of the constraints (§) and (0), we find that for Q = 1 

1 . . 
n f + n^ = -(1+t z ), 



n b + n 



x 



5 (l-r.). 



(13) 



For t z = 1 these equations revert to the constraints for a slave boson representation, when 
t z = — 1, they revert to those of a slave fermion representation, i.e an "up" superspin 
corresponds to a slave boson state, |(1 +r z )\ip) = \4>f)i a "down" superspin corresponds 
to a slave-fermion state |(1 — T z )\ip) = \iPb)- In the supersymmetric approach, a partition 
function of a Hamiltonian H, involves tracing over both slave boson and slave fermion 
representations, 



The trace over both subspaces means that the derived path integral has a U(l) x S77(l|l) 
symmetry and new dynamical degrees of freedom. In the slave fermion and slave boson 
schemes, Fermi liquid and magnetic phases are manifested as "Higgs phases" of the U(l) 
gauge group.0 The enlarged U(l) x SU(1\1) gauge group unifies the slave boson and slave 
fermion schemes, but also extends beyond it to furnish a potentially wider class of Higgs 
phases. For instance, suppose if is a Hamiltonian, such as the t — J model with both 
magnetic and paramagnetic phases, then we expect (r 2 ) = —1 in the anti-ferromagnetic 
(insulating) ground-state and (t z ) = +1 in the paramagnetic ground-state, but in addition, 
there is the possibility of new saddle-points where (t z ) lies between these two extreme values. 

We end with a discussion on the formulation of the t — J model as a supersymmetric 
large-iV expansion. To handle anti-ferromagnetic interactions and electron hopping in a 
large iV expansion, we adopt the Read-Sachdev scheme, using Hamiltonians that are globally 
invariant under the unitary symplectic group This group is a subgroup of SU(N) 

(defined only for even values of iV = 2n), so its generators are a subset of the Hubbard 
operators. Moreover, the groups SP(2) and SU(2) are equivalent. In SP(N), the spin 
components are divided into an equal number of "up" and "down" values a e (±1, . . .±iV/2); 
the unitary matrices of SP(N) satisfy the condition U T eU = e, where e aa ' = sgn(o")5 cr _ -/. 
The SP(N) t — J model is written^! 



Z= E Me 



s 



# = -^£[^o(*)X OCT a)+H.c] 



J 



i,3 3 



(14) 



where A/} = Z)o-^o-o-(j) is the number of particles. In the supersymmetric representation, 
this model becomes H + J2j Kj 



H = ~Tt YKfl<f>i + blxi)(4>}f j( r + xjM + H.c] 
iV (id) 

J 



(15) 



where Kj = Xj(Qj — Qo) + Cj(Xj ~ Yq) describes the constraints at site j, Mj = nf(j) + nb{j) 
and 



A; 



fiafja 1 fiabja' 

describes the singlet valence bonds between site % and site j. This Hamiltonian is invariant 
under the global SP(N) transformation and the local U(l) x S77(l|l) gauge group. The 
family of models with (Qo,Y ) = (N/2,0), (N even) are of particular interest. Two points 
deserve special mention: 

i) In a path integral treatment, by carrying out a local gauge transformation ipj — > gi(r)ipi 
and integrating over gj, one obtains a supersymmetric Lagrangianllll, C = C susy + H, where 



C 



susy 



JM 



d T - Xj + Qt 3 



This is the starting point for the study of the various Higgs phases of the model. In each of 
these phases, one of the fermi fields is absorbed into the fluctuations of the gauge field. For 
instance, in paramagnetic phases the slave boson condenses and by fixing 



V>j = 9 j 



ijo-N r i 



f . . . f '■ 

the slave fermions Xj are absorbed into the gauge field. Similary, the Schwinger boson field 
b a condenses in an ordered anti-ferromagnetic phase, absorbing a component of the f a fields. 



9 



More complex Higgs phases, in which fermi fields of the bond variables are absorbed into 
plaquet fermions also become possible. 

ii) The Lagrange multiplier Q which imposes the constraint on Yj gives rise to a self- 

It i 

consistently determined spin interaction Hi = — Q@j@j> resembling recent approaches to 
the Hubbard model in which spin interactions self-consistently renormalize to enforce local 
constraints^. The Gaussian fluctuations of the 9 fields associated with this spin interaction 
play a crucial role in enforcing the constraints between slave boson and slave fermion fields, 
and non-trivial results depend on the inclusion of these fluctuations in the effective action. 

In conclusion, we have presented a supersymmetric representation of Hubbard operators 
in which both the operators and the constraints are invariant under the action of the su- 
pergroup U{1) x S77(l|l). This approach avoids the need to choose between a fermionic, 
or bosonic representation for spins. The underlying U(l) x S77(l|l) gauge group is larger 
than the simple product of two U(l) gauge groups. Broken symmetry saddle points of 
this enlarged group provide the opportunity to study the interplay between magnetism and 
paramagnetism. 
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